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Abstract 

Quark mass effects are included in the calculation of polarized 
structure functions. In particular, the validity of fixed order pertur- 
bation theory and of massless evolution is studied in the framework of 
heavy quarks structure functions. The polarized version of the ACOT 
and MRRS interpolating schemes for the evolution of massive quarks 
distributions is also developed and studied. The different behaviours 
of the various approaches in x and are shown. 



1 Introduction 

There has been a great interest in ways to improve our knowledge of the 
heavy quark content of the structure functions, both theoretically ^ ^ 
^, ^, 1^ and experimentally [0]. It is an important research subject for two 
main reasons: 1) The gluon distribution in the proton has been extensively 
studied through heavy quark lepton production; 2) There are theoretical 
uncertainties on how to determine heavy quark distributions in the threshold 
region and how to resum the large logs containing the quark masses. 

So far, unpolarized electroproduction has been the major focus of atten- 
tion. The reason behind this activity is the large charm component of F2 
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measured at HERA However, future experiments dedicated to measure 
A^f P], the polarized gluon distribution of the proton, will also probe the 
polarized charm distribution through open charm production. Hence the im- 
portance of a proper treatment of massive quarks in polarized Deep Inelastic 
Scattering (DIS). In earlier attempts to estimate the charm component of 
Qip, attention was mainly restricted to its anomalous contribution |10|, |11 



and to the magnitude of Ac, the charm singlet axial charge |12, 13, 14 . In 



particular, in Ref. [|TT| the x dependence of the anomalous gluon contribu- 



tion to heavy quarks was calculated. One of the conclusions of that work 
was that the inclusion of charm was quite relevant to the magnitude of Ag. 

In any case, a detailed study of the x dependence and of the evolution 
of polarized heavy quarks is still missing. There has been an early work of 
Buza et al. where NLO corrections to the Photon Gluon Fussion (PGF) 
cross section is calculated and their x dependence studied. However, as it 
will be clear soon, it is desirable to expresses the heavy quark component of 
Qip in terms of heavy quark distributions. The main reason for this are the 
large logarithms which necessarily appear in the partonic cross sections. In 
this line, it is the purpose of this work to investigate how and when massive 
quarks can be treated as partons in polarized DIS. 

In massless QCD, the polarized structure function gip{x, Q"^) of the pro- 
ton can be written as: 

9ip{x,Q') = f-{C^%z.Q^^^)Aq^'{x/z,^^) 

J X Z 

1 -^i 

fc=l 

-rC,{z,Q\ii^)Ag{xlz,ii^)\}, (1) 
where the quark nonsinglet and singlet coefficient functions are 



while the gluon coefficient function is 



'^(z, Q^li?) = 5{z - 1) + Q^/^^) + (2) 



C,{z,Q'/^?) = ^C«(.,QV/i^) + (3) 
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and the sum runs over the hght flavours. The singlet quark distribution is 



A^{z,^i')=J2Aq,{z,f,'), (4) 

i=l 

and the nonsinglet quark distribution is 

^QNS = i E fe,^ - - E el] Ag,(z, /i^). (5) 



We set the renormalization scale to be equal to the factorization scale, /j,^. 
The dependence appearing in Eqs. (|^) and (||) is in the form of a 

logarithm. Given the x dependence of the axial charges in Eqs. (§) and (|^) 
at some scale yU^, we can calculate the structure function at any other scale 
with the help of the perturbatively calculated coefficient functions and 
splitting functions. Basically, there are two ways to perform this evolution. 
First, we could fix fi"^ and do the evolution via the logarithms appearing in the 
coefficient functions, a procedure which is called Fixed Order Perturbation 
Theory (FOPT) pO|, |21[. This method becomes problematic when » fi^ 



because large logarithms appear. Indeed, Cg in Eq. (|]) is proportional to 
as(/i^)/n((l — z)Q'^/zfi'^) and the LO approximation will be reliable as long 
as this term is much smaller than unity. On the other hand, in FOPT the 
structure functions are determined to any order once the relevant partonic 
cross section is calculated to that order. It is not necessary to know all the 
splitting functions either to evolve the parton distributions. Otherwise, we 
could set /i^ to the natural scale and do the evolution via the DGLAP 
equations for Aq^^{z,Q^), AT,{z,Q^), and Ag{z,Q'^). In other words, the 
logarithms in would be resummed. This procedure would be the most 
adequate as in the case of FOPT we have always the problem of not knowing 
up to which the evolution through the logarithms can be trusted. 

The situation is more involved when we deal with massive quarks. To 
begin with, we have an extra scale in the coefficient functions, the heavy 
quark mass m/j. For >> ml, the mass sole effect would be to replace 
yU^ by ml in the logarithms of the coefficient functions, which means that 
the quark mass is being used as the infrared regulator. The corresponding 
large logarithms can then be resummed in a similar fashion to the massless 
quarks. Indeed, this has been done to all orders of perturbation theory by 
Buza et al. B for the unpolarized structure functions. When ~ m^. 
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this resummation can not be done. It then seems natural to use FOPT to 
calculate the structure function in the threshold region, using the full mass 
dependent form for the coefficient function. However, it is desirable to have 
a scheme that interpolates the two regions. The recently proposed methods 
of Aivazis et al. (ACOT) and of Martin et al. (MRRS) [g do this inter- 
polation. In the MRRS scheme, the relevant splitting functions are modified 
in the threshold region in order to incorporate mass effects in the evolution 
of the massive quark distributions. It means that a heavy quark distribu- 
tion is introduced at threshold which is evolved with the help of the relevant 
mass modified splitting functions. Moreover, the gluon coefficient function is 
modified to avoid double counting. When » m|, their expressions are 
reduced to the massless convolution formulas. In the ACOT formalism, the 
splitting functions preserve their massless form but the coefficient functions 
retain their full dependence in m\ to all orders of perturbation theory. Inci- 
dentally, we will see in the present work that the ACOT scheme is natural 
in polarized DIS. 

It is the aim of this paper to study the evolution of the heavy quark com- 
ponent of gip using three different approaches. In section 2, the basics of the 
FOPT and of the massless evolution of gip{x, Q'^) are spelled out. In section 
3 we calculate the modifications that the quark masses induce in the polar- 
ized splitting and coefficient functions. We then implement the interpolation 
schemes proposed in |]. The different methods for calculating gip{x,Q^) 
are compared in section 4, and their validity is discussed in section 5. 

2 g^p{x, Q^) in FOPT and in Massless Evolu- 
tion 

In the absence of intrinsic heavy quarks in the proton, the dominant pro- 
cess for heavy quark production is PGF. This partonic cross section starts 
at order a^, and once it is known, one can in principle calculate the corre- 
sponding structure function at any scale. Hence, if at the scale /x^ the heavy 
quark distribution is zero, then the heavy quark component of the polarized 
structure function gip{x, Q"^) to order ag is given by: 

9\,M')=el^^ f'-H^'\Q\mlz)Agi^,^x/z), (6) 

47r J ax Z 
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where the PGF cross section for a massive particle, H^^\ is |I5| , |I6[: 



H('\Q\ ml z) = Tf[A{2z - l)ln (^i±^^ + 4(3 - Az)sq], (7) 
with Tf = 1/2 and 



sq 



Ami 



\ (1 - z)Q^ 



(8) 



The factor a = 1 + Aml/Q'^ in (|) ensures that the photon carries enough 
energy for the resolution of a heavy quark pair. 

The questions to be answered are: given the scale /i^, how good is Eq. 
(H) to describe gip{x,Q^)7 Up to which x and is Eq. valid? In the 
absence of data, the answer to these questions should be in the comparison 
with different theoretical approaches. The main constraint to Eq. (j^) is 
as{^'^)l'n{Q'^ /ml) << 1, a condition that is fulfilled in the threshold region 
(Q^ ~ ml). However, when Q'^ » ml Eq. is reduced to: 

i/«(g^ ml z) = Tj[A{2z - l)ln f + 4(3 - 4z)], (9) 

and the large logarithms, mentioned in the introduction, appear. In that 
case, FOPT is not reliable any more. In an early study ||2^ it was argued 



that Q"^ /ml < 100 would be a good approximation to FOPT. As we will 
see later, this is true only in the x region studied in that work (x > 0.01). 
The discrepancy between FOPT and DGLAP evolution in fact grows as x 
becomes smaller, becoming as different as 100% already at x = 10"'^. 

Of course, one could circumvent the problem with the large logarithms 
assuming massless evolution for heavy quarks over the entire region. This 
approximation has been used in the literature for the unpolarized heavy quark 



distributions [|17|, 0. In such heavy quark distribution function is 

introduced at a given scale, say ml. Its value at any other scale is then 
predicted by the usual DGLAP equations. In this approach, the heavy quarks 
are resummed as massless particles and Qi would be given by: 
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J X Z 

+ ^C7«(.)A,(a:/.,Q^)], 



(10) 



with Ah{z, ml) = and C^^^ given by Eq. (^ taking = ml. The massless 
quark coefficient at 0{as) is 



with m, 20 



C,(z,a(Q^))=5(z-l) + ^^C«(.) 



(11) 



l + z^ 
I- z 



'ln{l-zY 
1 - z 



1 - z 



2{l + z)ln{l ~ z) 



Inz + 4 + 2^ + 5(1 - 2)(-4C(2) - 9) 



(12) 



Although massless evolution would give the correct answer at >> Tnl, 
it does not solve the problem of evolution of heavy quark distributions. To 
treat this problem, we apply the interpolating method of Refs. [|T], ^ to the 
polarized case. 



3 in an Interpolating Scheme 

As massless evolution for massive particles would be clearly wrong when 
~ ml, and FOPT would present problems at very large values of Q^, it is 
desirable to have an scheme that interpolates between both. This procedure 
should reproduce Eq. (H) at the threshold and Eq. (p!OD when » ml- 
Despite the fact that the ACOT scheme |jl| does this interpolation, in a 
recent work Martin et al. |^ developed an alternative scheme where such 
interpolation is also obtained. The main ingredients of their work are: 1) 
A heavy quark distribution is introduced at threshold, and its evolution is 
according to DGLAP type equations, with mass modified splitting functions; 
2) The PGF cross section is modified to account for similar contributions 
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coming from the evolution of the heavy quark distribution. This treatment 
ensures the absence of double counting, in the same spirit of the ACOT work. 

As in our work we are interested in the polarized heavy quark component 
of the structure functions, we have to consider the mass modifications in the 
polarized splitting and coefficient functions in order to implement the MRRS 
scheme. We start computing the polarized gluon splitting function APhg in 
LO. The decay probability of the gluon into a polarized quark pair, where 
the quark carries a fraction z of the gluon momentum is given by: 

dPhg{z)dz = -^APhgdzdln{pl + ml) (13) 

ZTT 

with the polarized splitting function being defined as: 

A „ / N — Z) A , 

^P.M = (") 

As usual, the four vectors will be defined in the Infinite Momentum Frame 
(IMF). The amplitude A for the decay, where the parent gluon is taken to 
be always in a state of positive helicity, is: 

A±+ = Tr[{h + m,)7^(/^^ - mnh,^^]erel. (15) 
The gluon polarization vector is: 

e± = (0,0, 71/2, ±^^), (16) 
and the quark vectors are 

h^(zp + ^^^,zp,p^), (17) 
2zp 

with a similar expression for k-j^, just z replaced hy 1 — z. We also use the 
convention e^^"^^ — +1. After a straightforward calculation, we get: 

p _ ^ ^(1 - ^) ^++ 

^h+g+ - J-f 



2 {pl + mir 
Tf(z' + -T^zil - z)] , 



Pn-g^ = T,[{l-zr + ^f^^z{l-z)j. (18) 



7 



The splitting function is then: 

AP,, = P,+,+ - P,_,+ = ]^{z' - (1 - zf), (19) 

which is protected against mass effects. That means that a decay of a gluon 
into a massive quark pair does not produce any extra breaking of hehcity. 
We can also get the unpolarized splitting function: 



1 Tn 
Pn, = Ph-,,+ + Ph-,+ = -iz' + (1 - zf) + ^-^zil - z), (20) 

z P±-r 

which agrees with the expression calculated in |^. 

Our next step is to calculate the APhh splitting function: In this case, 

AP,, = c/^\ '^i + mi^(l-^)2' (21) 



with 



A±+ = Trlih + m,)^4^7.(^;. + m,)i4^7M]e* • e. (22) 



The scalar product of the gluon polarization vectors is 6^^ — kgh^lkgl"^, with 
i,j = 1,2,3., and kg the gluon four momentum vector in the IMF. In the 
calculation, we took the initial quark to be always in a state of positive 
helicity. For the helicity conserving part, we get exactly the same 

result of the massless (and also unpolarized) case. For the helicity violating 
amplitude, we get y4_+ = — 4m|. It results in: 



W = -2C.--(l-.)^^^^|i-^ (23) 



Hence, 



AP 



hh 



.1-2 

A similar calculation shows that 



l-z + -S{l-z)+2z{l-z) 



mt 



2''~ ~' ' ">i + m^(l-z)2 



(24) 
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mi 



Pg+n+ = Cf[--z{1- z) \ ^ , (25) 



and 



P,-H^ = Cf (^-^—^ - z{l - z) , f\ \ . (26) 

Thus l^Pgh is also protected from mass effects. The same is not true for its 
unpolarized counterpart. 

The last quantity to be calculated is the polarized partonic cross section, 
gf^ . It can be extracted from the Born approximation to the antisymmetric 
part of the partonic tensor. 

The general structure for the antisymmetric part of the partonic tensor 

is: 

0^7 = -^e^-'-'qalsM^, Q') + (s, - '-^] g^iz, Q% (27) 
2p.q \ p.qj 

where p is the parton four momentum vector, s its spin vector and q the 
photon four momentum vector. If we define a general projector P^'^ = 
e'^'^'^AaBfj, such that: 



P'-'^i = -^{A.qB.s - A.sB.q)gf\ (28) 
^ zp.q 



then 



9''' = , ^ \ ^ P^'-i- (29) 

nih A.qB.s — A.sB.q ^ 

In the Born approximation, the antisymmetric part of the partonic tensor is 



A 



m„£^.a/3sV5(2p.g + g'). (30) 



The use of Eq. (pO]) in Eq. (pQ]) then renders a coefficient function in LO 
which is unmodified by the presence of quark masses: 

# = 5(^-1). (31) 

We see that in LO, the polarized splitting functions and the coefficient 
function are protected against mass corrections, with the exception of the 
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AP/i/j term. However, as shown in 0, any mass effect appearing in the LO 
expressions, will be felt only in NLO. This can be seen from the following. 
The general expression for the decay of a parton in LO would be: 



dP{z)dz = \ + 9{z)j:t—-^, dzdpl. (32) 



«M1 [ /(^) + ^(^) af. 

27r \ pi + fn"^ (pj_ + fn^ 

As is well known, if we integrate Eq. (^) n times in p']_ using = 0, with a 
strong ordering in p^, and sum over n, we obtain the [a^fi"^) / a{Q^)]^ of the 
renormalization group equations, where P is the splitting function. When 

7^ 0, logarithms in a are systematically lost, which is a typical feature of 
NLO resummation. Besides this, the integrals of the mass terms in Eq. (p^) 
are problematic, non analytical, and resummation is impossible. In the same 
way, any mass modification of the NLO splitting functions, would appear 
only in a NNLO study. With this in mind, we can then make a fully NLO 
evolution of a polarized heavy quark distribution using the same massless 
expressions, the only exception being the LO APhh term. 

The next step is to define the massive PGF cross section which interpo- 
lates the massless coefficient, Cg in Eq. (^), and the full massive coefficient, 
Hg in Eq. (P). As noticed in Ref. a part of Hg has to be subtracted 
because when evolving the heavy quark distribution, the APhg splitting func- 
tion will also create heavy pairs. In the massless evolution such a subtraction 
is already implicit in the absence of the ln{Q'^) term. Using the slow evolu- 
tion variable, ^ = (1 + m\/Q'^)x = bx, as defined by Brock ||22|, we write the 
heavy quark component as: 

Ah{x) (g) g[°\x) = ^5iz - l)Ah 0^ = Ah{0. (33) 

Integrating the gluon splitting part of the LO DGLAP equation for Ah, we 
obtain: 



47r Jf z JQ^ . 2 

> ^ mm 

If we use asAg as a scale invariant, then Eq. (|3^ ) is reduced to: 

Aft«) = ^ta(5)£fA<,(.,Q^)^^f£M. (35) 
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where we use m| for the minimum value of the scale. According to the APhg 
splitting function, we should start to evolve the heavy quark distribution at 
m|. This is explicit in Eqs. (|T^) and (|55D. And we actually implement this 
scale in our calculations. However, to avoid any unconsistency with the scale 
necessary for heavy quark to be resolved as a physical object (being probed 
by the photon), we will always calculate g^^ at scales where > ^rn\. 
We then incorporate into the heavy quark coefficient function a threshold 
factor of the form 6{1 — 4m^/Q^). That means that below 4m^, gip{x, Q"^) is 
given by PGF only. As the choice of the threshold factor is not unique, the 
calculation of the heavy quark structure functions is, in a sense, ill defined. 
We return to this point in the next section. 

Finally, we write down the convolution formula for the heavy quark com- 
ponent of the polarized structure function: 

9l,{^,Q') = el ['-C,{z,Q')Ah(^-^,QA 

Jbx Z \ Z J 

+ el^^ f'^HW(z,mlQ^)AgC-,Q^) 

47r Jax Z Z 

_ (36) 

47r ym^ ) Jbx z 2 \ z ) 

When we take » ml the last two lines of Eq. ( p6| ) are reduced to the 
expression for massless quark production through PGF, and we get Eq. (p!o|) 
corresponding to the massless case. We notice that this is an exact result. 

In light of the resulting Eq. (^6]), it is instructive to analize the short- 
comings of both MRRS and ACOT approaches for the calculation of the 
heavy quark structure functions. The most evident and serious problem 
with MRRS is the fact that one does not recover the massless, MS calcu- 
lated, unpolarized gluon coefficient function when the limit » m| is 
taken 0. Hence, it can be hardly claimed that the defined heavy quark 
distributions are MS distributions. In the same line, MS'-like schemes are 
mass independent, meaning that evolution can not possibly depend on the 
mass of the quark, a requirement not fulfilled by the MRRS prescription. 
However, because of the mass independence of both the polarized quark co- 
efficient function and the polarized splitting functions in LO (the exception 
being APhh but it has no numerical relevance - see next section), Eq. ( p^ ) 
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is really the ACOT version for the calculation of the polarized heavy quark 
structure function. Hence, opposite to the unpolarized case, the polarized 
gluon coefficient function has the correct MS form in the limit » m\. 
In a sense, it can be said that the ACOT formulation is preferable over the 
MRRS formulation for gip{x,Q'^). In any case, the ACOT formulation also 
has formal drawbacks. Specifically, the gluonic contribution to the structure 
function, which is of order a^, enters together with the a° quark contribution 
only, when the as quark corrections should also be considered. However, as 
stressed in the ACOT work fll], the as quark corrections are numerically in- 
significant when compared to the gluon contributions. In the particular case 
of the present work, the ag quark corrections were estimated, and their con- 
tributions are negligible. Finally, both schemes have ambiguities related to 
the choice of the factorization scale /x^, when physical observables, like struc- 
ture functions, should not depend on the chosen value for yU^. Fortunately, 
in practical calculations this dependence is not a real concern 0. 



4 Results 

We are now ready to calculate gip{x, Q^). Although we will restrict ourselves 
to charm, the extension to other heavy quarks is immediate. We will use 



the NLO GRSV [^] "standard scenario" parametrizations for the gluon and 
quark singlet distributions as the input in our calculation. We made this 
choice because of the low scale where this parametrization is supposed to be 
valid, which allows us to evolve through charm threshold. The strategy used 
for evolution was the following. We first evolved the gluon and the singlet 
quark distributions from 0.34 GeV^ to m^, with 3 active fiavours. From ml to 
Q^, the charm quark distribution is obtained from the evolution of the singlet 
quark distribution with 4 fiavours, using as input distributions AE(2;, ml) = 
and the Ag{z,ml) obtained from the evolution from 0.34 GeV^. The charm 
distribution at is then Ac{z,Q^) = AT,{z,Q^)/4:. The gluon distribution 
at is obtained from the singlet evolution of the gluon distribution with 
four active fiavours, using as input distributions AIl{z,ml) and Ag{z,ml), 
which were obtained before when evolving from 0.34 GeV^. Once we obtain 
the gluon and charm distributions at Q^, we can use the various convolution 
formulas, Eq. (P) for the pure PGF evolution via FOPT, Eq. (|36D for the 
improved MRRS method, and Eq. (^) for massless evolution. 
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We use rric = 1.5 GeV for the charm mass, and ^qcd — 0-2 GeV in order 
to agree with GRSV. The couphng constant at any scale, from ij? to Q^, is 
given by the solution of the NLO equation for a^: 

asm ^ 1 _ A MMQVA|c^ 

We notice that the use of this equation at low values of such as 0.34 GeV'^ 
can lead to difficulties. In this case, the use of the transcendental equation 
for as from which Eq. (|37D is derived, would be more appropriate. However, 
to be consistent with the way the parametrizations were built, Eq. ( |T7D have 
to be used. 

In Fig. ID we present the axial singlet charm distribution of the proton 
at 10 GeV^. It is obtained as explained before, evolving the singlet distribu- 
tion from m^. This result shown is independent of the mass modified ^Pqq 
splitting function. It means that even if AP^^ is the only term which receives 
any mass modification, it is also protected against mass terms for practical 
purposes. It is possible to fit the charm distribution as: 



xAc(x, 10 GeK^) = x°-^^(-0.055(l-x)^ + 0.206(1 -x)^) 

+ x°-^^(-0.286(l + 0.261(1 -x)^). (38) 



The particular form for Eq. (^) was inspired by Brodsky et al. ||2J] consid- 
erations on the behaviour of the polarized distributions. The integral over x 
of Eq. (ID is 

Ac(10 GeV"^) ~ 2.6 x 10-^ (39) 

which is completely irrelevant. However, we notice that charm is allowed in 
a large portion of the existing data for gip, meaning that the study of the x 
dependence is relevant, mainly in the small x region. 

In Fig. Hthe charm component of ^^ip at 10 GeV^ is presented. The dotted 
line is the result from FOPT, through Eq. (^. In this 2.25 GeV^ 

and as is calculated using rij = 4. The dashed line is the result from massless 
evolution. The continuous line is obtained with the help of Eq. (^), which 
is the scheme where low and high are consistently incorporated. We 
see that down to a; ~ 4 x 10~^, the three methods of calculating g^Jx, Q^) 
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Figure 1: The x dependence of the charm singlet axial charge at 10 GeV^ . 

give very close results. However, as x gets smaller FOPT fails to follow the 
other two methods. We think this failure happens because at x = 10""^, 
as(/i^)/n((l — x)Q'^ /xml) ~ 3, and the NLO corrections calculated in 
would need to be incorporated. 

On the other hand, we plot in Fig. ^ the x dependence of the pure PGF 
term, second line of Eq. ( ^61) and the pure charm quark distribution, the first 
plus the third lines in Eq. (|36|). The dashed line is the PGF contribution 
and the continuous line the charm contribution. Comparing with Fig. ^ 
we see that the PGF alone would reproduce the FOPT calculation, which 
is a natural result as Eq. (^) is pure PGF. But the charm contribution 
is almost as large as the PGF, hence the discrepancy between the dotted 
and continuous hne of Fig. ^ That happens because, as stated in the 
end of section 3, we choose not to add any ad hoc threshold factor in the 
charm coefficient functions. If we had added, for instance, the factor / = 



X 
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-5 ^ i 

0.0001 0.0010 0.0100 0.1000 

X 



Figure 2: The charm component of gip{x, Q"^) at 10 GeV^. The dotted hne is 
the result from FOPT. The dashed hne is the resuh from massless evolution 
and the continuous hne represents Eq. (|36|). 

(1 — 4:ml/Q'^)6{l — Aml/Q"^) as it is done in 0, the charm quark component 
of gip{x, Q"^) shown in Fig. ^ would be dramatically reduced. In this case, the 
line of Fig. |^ would be very close to the FOPT result. However, the factor / is 
arbitrary and other choices for it would produce different contributions from 
the charm evolution. In this sense, the calculation of glp{x, Q^) at threshold 
is ambiguous. In our case, we made the choice that if charm can be resolved, 
then it is given directly by Eq. (p6|). In other words, our threshold factor is 
the 6 function only. 

We show in Fig. ^the result of evolving the singlet part of gip{x, Q^) with 
three and four flavours. The procedure used was as the same used before. For 
three flavours, we evolve directly from 0.34 GeV^ to 10 GeV^, using n/ = 3 
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0.0001 0.0010 0.0100 0.1000 
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Figure 3: The continuous line is the charm component of Eq. while the 
dashed line is the PGF term only. Both curves are calculated at 10 GeV^. 

(s) 

and ^QCD — 0.248 GeV . For the four flavours component, we evolve from 
0.34 GeV^ to ml in the three flavour theory and from ml to 10 GeV^ in the 
four flavour theory. Of course, the procedure is not completely consistent 
as the GRSV parametrization was built using rij = 3 only. However, it is 
expected that at the parametrization is valid, meaning that the results 
shown here really point to the modifications that a fourth flavour would 
induce to gf^i^x^Q"^). Indeed, we see that at x = 10~^, a difference of the 
order of 30% is present. 

Finally, Fig. ^ shows glp{x, Q^) at 100 GeV^. As expected, at this scale 
the improved equation for massive convolution reproduces exactly the mass- 
less approach. FOPT, however, presents the same sort of behaviour shown 
at 10 GeV^. The reason for this behaviour is the same as explained before. 
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Figure 4: The singlet part of gip{x,Q'^). The continuous curve was evolved 
with 4 flavours, while the dashed curve was evolved with 3 flavours. Both 
curves were calculated at 10 GeV^. 



5 Final Remarks and Discussion 

We have computed the charm component of the polarized structure function 
of the proton, gip{x,Q'^). We used three different methods in these calcu- 
lations and it is remarkable that already at 10 GeV^, massless evolution is 
a very good approximation. Of course, this conclusion is blurred by the 
umbiguous appearance of a possible threshold / factor in the heavy quark 
coefficient functions. In any case, it is seen that FOPT is close to massless 
evolution down to a; = 10~^, at least up to = 100 GeV^. This conclusion 



complements an early observation where it was argued that FOPT and 
DGLAP evolution were equally good. In truth, we see that when x gets 
smaller, the large logarithms associated with the scale are enhanced by the 
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Figure 5: The same convention of Fig. 2, calculated now at 100 GeV"^. 

1/x factors, something which does not happen in the usual evolution ap- 
proach or in the improved method of Eq. (|3^) because the scale logarithms 
are absent. FOPT then breaks down. Moreover, we saw in Fig. ^ how impor- 
tant the inclusion of a heavy flavour can be for the singlet component of gip, 
even if the amount of Ac in the proton is small. Finally, when calculating the 
mass corrections to the splitting functions, we saw that they are protected 
against mass effects. Thus, effectively, Eq. ( ^6|) is the polarized version of 
the ACOT scheme ||l| for heavy quarks. 
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